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Introduced recently, the concept of hierarchical degree allows a more complete characterization of 
the topological context of a node in a complex network than the traditional node degree. This article 
presents analytical characterization and studies of the density of hierarchical degrees in random 
and scale free networks. The obtained results allowed the identification of a hierarchy-dependent 
power law for the degrees of nodes in random complex networks, with Poisson density for the first 
hierarchical degree (obtained through master equation approach). Exact results were obtained for 
the second hierarchical degree in scale free networks. 

PACS numbers: 89.75.Fb, 89. 75. Da, 02.10.Ox 



I. INTRODUCTION 

The characterization and analysis of complex net- 
works [l|, |2j involves the use of measurements capable of 
expressing important properties of the network. While 
traditional features such as the node degree and cluster- 
ing coefficient have been successfully used for that pur- 
pose, such measurements are intrinsically local, in the 
sense that they are affected only by the immediate neigh- 
borhood of each reference node. Among more global mea- 
surements, such as the network diameter and the shortest 
path between any two network nodes, the concept of hi- 
erarchical node degree has been introduced recently and 
considered for the characterization and analysis of the 
connectivity of complex networks 0, IE IS ■ Given a 
reference node i, it is possible to define the k — th relative 
hierarchical level as corresponding to the network nodes 
which are exactly at minimal distance of k edges from 
node i. The hierarchical node degree at level k corre- 
sponds to the number of edges between the hierarchical 
levels k and fc+1, therefore defining a natural hierarchical 
extension of the concept of node degree. Note that the 
hierarchical degree, starting at level 1 (which includes the 
immediate neighbors of i) tends first to increase with the 
hierarchical level k and then to decrease as a consequence 
of the finite size of the network Q ■ Because the hierar- 
chical node degree is a global measurement (i.e. it con- 
siders, in hierarchical fashion, all nodes in the network), 
it provides more information about the connectivity of 
the reference node with the remainder of the network. 
For instance, the number of hierarchical levels which are 
required in order to contain a percentage (e.g. 50%) of 
the network nodes, provides a valuable information about 
the accessibility of that reference node. The smaller this 
number of levels, the more connected the reference node 
is with the rest of the network. Several other topologi- 
cal properties characterizing the context of the reference 
node can be identified by using the hierarchical node de- 
gree and other related features 00]. 



An interesting question implied by the introduction of 
the concept of hierarchical node degree regards the an- 
alytical characterization of such features expected from 
random and scale free networks. The current work ad- 
dresses such an issue by using mean field and the master 
equation methodologies E| • It is shown that a power 
law is obtained for the hierarchical degrees of nodes in 
random networks, with Poisson density being obtained 
for the first hierarchical degree. We verified that the 
second order hierarchical degree for scale free networks 
involves a logarithmic correction of the first order degree 
studied in 0, 0] . In addition, we obtain exact results 
for the second order hierarchical degree distributions in 
scale free networks. 



A. Hierarchy and hierarchical degree 



A set of hierarchy h is defined as that containing all the 
vertices which are located at a distance of h connections 
from a reference vertex. The hierarchical degree, or su- 
perior order degree, is the measurement of the number of 
connections in a set with a determined value of hierarchy. 
For ft, = 1, we have the first order degree which is well 
known in the literature 0, • I n this article we are 
intended to study the behavior of graphs where h > 2. 
We will denote fc z as the degree of order h of the refer- 
ence vertex i; and the set of hierarchy h of the reference 

vertex i will be denoted as {i[ h \ 4 , 4 > •••> Mfc-i) }■ We 

will also denote as the number of vertices having 

hierachical degree h equal to x, and 2V a ,&,—,n wm " denote 
the number of vertices presenting the first order degree 
as a, the second order degree as b and so on. 
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II. DEGREE EVOLUTION OF A VERTEX 

Mean field theory is widely applied in complex net- 
works research 0, 0, . Through this boarding, the de- 
gree evolution of the particular vertice is proportional its 
connection probability (kernel function - Pj(t)). If this 
probability has a complex form, the equations become 
non-trivial. 

Generally, the hierarchical degree h depends on the 
probability of connection of the vertices belonging to this 
set of hierarchy. Thus, we can write: 



dkf\t) 



u(h-l) 



(1) 



In the following we consider random graphs where 
Pj(t) is the same for all vertices, as well as scale-free 
networks proposed by Barabasi- Albert 1] where Pj (t) oc 



III. DISTRIBUTION OF THE SECOND ORDER 
CONNECTIVITY DEGREE 

In this section we will study the connectivity degree 
distribution, which is an intrinsic characteristic of the 
complex networks and is capable of distinguishing all 
types of graphs. Henceforth, for simplicity's sake, the 
probability of vertex j to make a connection will be rep- 
resented as V. (i) (t) and not Pj(t). 



A. Initial States and Structures 

A vertex with first order degree fcW and second or- 
der degree k^ defines one or more structures which are 
topologically identic and are denoted by (k^\ k^). For 
example, Figure Q shows all possible structures (3,3). 

An structure has one or more initial states. An struc- 
ture (x\,X2) is the initial state of another structure 
(jji, 2/2) if the following condition is met: 



(xi +Si,x 2 +S 2 ) —> (2/1,2/2) 



where: 



Si 



1 



(2) 



(3) 



Therefore, two processes exist which take one structure 
into the other: 



(.xx + 1,2/2) (2/1,2/2) (4) 
(2/1,2:2 +1) -> (j/1,2/2) (5) 

The initial states of structure (3, 3), given by (2, 3) and 
(3, 2) are represented in Figure |3 



B. Describing the Processes through the Master 
Equation 

The probability of transition of the initial states in 
a structure give us the production rate of that struc- 
ture. For example, let us study the formation of struc- 
tures (fcW, fc( 2 )). As we know, their initial states are: 
(fcW - l,kW) and (k^,^ - 1). Then, the transition 
probabilities are given by: 



^ ) (fc( 1 )-i,fe( 2 ))->(fe( 1 ),fe( 2 )) — PkW-l 



^ > (fc( 1 ),fc(2)_i)_ > (fc(i) i fc(2)) — Vi-a 



Therefore, the formation rate of the structure 



are: 



N 



— 7 : ' fe (i)_ 1 A fc (i)_ lfc (2) 
+ S lk (i)V k {2)N^ {2) 



k (l) +k (2)_ 1 N k (l) k (2)_ 1 



The transition probabilities of structure (k^\ k^) is: 

P{ fc(i), km)^(kW+iM 2 '>) = PkW 

P{ fed), fc( 2 ))->(fc(i),fc( 2 )+i) = ^(i)+fe< 2 ) 

Therefore, the vanishing rate of the structure 
(fcW,fc( 2 )) is: 

N_ k (l) k (2) = [V k (l) + V k (l) +k {2)]N k (l) k (2) 

The master equation governing the evolution of 

N k (i )k (2) is as follows: 



dN k g) k (2) 
dt 



- [PkW +Pk( 1 )+k( 2 )] N k( 1 )k< 2 ) 

+ 5 lkW V kW N£l (6) 



In the following we will show the solution of the equa- 
tion above for random graphs as well as for scale-free 
networks. 



IV. APPLICATIONS 

A. Random Graphs 

In this work we use non-static version of non-directed 
random graphs. Instead of generating a network from a 
symmetric random matrix, we define the following pro- 
cedure: Beginning with Af disconnected points, we select 
two vertices in the network at a time and, with proba- 
bility p, we establish a connection between them. Con- 
nections linking a vertex to itself are not allowed, that is 
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FIG. 2: Initial states of structure (3,3). 



to say, the two selected vertices must be different. The 
probability of a vertex j to be selected at a time t is given 
by: 



As is known, N^\t = 0) = Af and Nj{t = 0) = 0,Vj. 
Solving this equation for — 0, we have: 



A> Af Af —1 N{N-lf 



(7) 



N^(t) = Afe 
And for k^ 1 ' = 1, we have: 



1. Degree Evolution of a Vertex 

Since Pj(t) is constant, we will denote it as it. There- 
fore: 



7T = P j (t) (8) 

By substituting Eq.(8) in Eq.(l), we have: 



dkf\t) 



(9) 



Assuming that kS — ght h , we find gh = n h /h\ and 



(10) 



2. First Order Degree Distribution 

The variation of the number of vertices with first order 
degree given as fcW can be described through the master 
equation: 



^ 7r L iV fe(i)-i Jv fc(i)J 



iV^^t) = nATte-*' = 7rtJV w (t) 

For other values of fc^ 1 - 1 , we conclude that the solu- 
tion of the master equation obeys the following formation 
rule: 



The solution for this recursion is: 



feW! 

The distribution of the first order degree, given by 



(Trt) 



feWg-TT* 



(11) 



This is the Poisson distribution with mean value 7rf. 
Figure |3| shows the time evolution of the number of ver- 
tices with k^ — 0, 1, 2 in a random graph with AT = 100 
and p — 0.9. 
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FIG. 3: Evolution of the number of vertices with first order degree fe' 1 ' = {0, 1,2} in a random graph with J\f = 100 and 
p = 0.9. 



B. Scale-free Networks 

The concept of scale free networks has been introduced 
by Barabasi and Albert and refers to graphs whose first 
order degree distribution follows a power law Q: N k m = 
A^ 1 )] - ". Most scale free networks found in nature have 
v in the interval between 2 and 3. The probability of 
connection in this model is given by: 



„(i) 



(i) 



= ir- (12) 

2t 



Also: 



k (h-i) + k (h) k (h-i) + k (h) 
V k w{t) = — t = ^ (13) 



2t 



1. Evolution of Vertex Degree 
Substituting Eq.(12) in Eq.(l), we have: 



dt 



2t 



(i) 

3 

2t 



(14) 



;From the definition of hierarchical degree, we have: 



,<h-i) 



,f) +k f-v = 8lhk f) + j2 q 



(i) 



(15) 



Therefore: 



dkf\t) 



dt 2t\ 3 3 



(h-i) 



(16) 



For h = 1 and the following condition: k^(t = tj) = 1 
we obtain the solution given by Barabasi- Albert in |3j: 

(17) 



*i l) ('> = i/r 

For h = 2, we use the integrating factor i -3 / 2 and we 
obtain: 



kf\t) = [ft+lnVt] kf\t) 

If we consider that, in instant where the new vertice 
is enclosed to network its average second order degree 

(2) (2) (2) 

is given by Kj , then: kj (t = tj) = Kj . From this 
condition, we have: 



kf\t) 



■5 
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FIG. 4: Values of kS in terms of j obtained through numerical simulations are adjusted by the curve k\ ' = 1 + hi j. 
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FIG. 5: Evolution of kgg (t) for a random graph with size Af = 100 and connection probability p = 0.9 and Scale-Free network. 
The curves intercept at t = r, indicating that at this growth stage vertex j in the scale free network has the same second order 
degree as vertex j in the random network. The log-log plot in the inset shows how r(tj) depends of tj, defining power laws 

T{tj)0CtJ a . 



Through numerical simulations, we can determine the 
value of n^p in terms of j. From Fig(4) we note that 

(2) 

Kj = 1 4 In j. Thus, the solution for the second order 
degree variation is given by: 



(t) 



(18) 



Eq.(15) is not simple to be solved but, using the same 
integrating factor as above, we can obtain a solution as 
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an integral recursive equation, as follows: 



kf )(t) = p h kf\t) + \Vt f kf-^ridt (19) 



2. First Order Degree Distribution 

In 01 j Krapvisky establishes the following master 
equation which governs the evolution of N^ l} , that is, 
the average number of vertices with first order degree 
feW: 



(i) 



<9< 



M(t) 



i.(i)AT (1) 



+ <5 



lfcCD 

The first term takes into account the probability that 
the new vertex connects to one of the first order degree 
equal to k^- 1 ' — 1. The second term expresses the probabil- 
ity of connection between a vertex with fcW. The delta 
function accounts for the fact that every added vertex 
has degree 1. By making the normalization M(t) = 2t, 
the solution of Eq.(15) is given as: 



k(1>[> + l)(fe(i) +2) 



it 



(21) 



The distribution of the first order degree, i.e. 
P S (*) = N S) (*)/^T, is given as: 



this quantity involves the solution of Eq.(6). Assuming 
a linear solution of the type: N k (i )k (2) — n k (i )k (2)t, we 
have: 



n fc(i)fc(2) 



7 ? fc ci)_irt fc (i) _ lfc (2) + 7 ? fc (i) +fc (2)_ 1 n fc (i) fc (2). 



Pi + Pl+kW + 1/t 
Substituting Eq.(13) into above expression: 



n fe(l)fe(2) 



(fc^ - l)n fc (i)_ lfc( 2) 
2 fc(i) + fc(2) + 2 

jfe(l) + fc(2) _ 1 



/ jfc(l) + fc(2) _ 1 \ 

+ Ufc (1) +^ (2) +2j nfe<1)fc(2) - 1 
'fc (2) n fc < 2 > ' 



4 + fc(2) 

For simplicity's sake, we make 



A u — 
F u = 



' (fc ( "-iKw-h 

2fc« +u + 2 

fcW +u-l 
2fcW + M + 2 



(23) 
(24) 



4^ J to (25) 



(26) 



/^(fct 1 ) +l)(fc« +2) 



(22) 



Figure |3] shows the evolution of k$ (t) for a random 
graph with N = 100 nodes and connection probability 
p = 0.9 as well as the same measurement for a scale free 
network. An intersection point can be observed at t = r, 
meaning that at this growth stage the vertex j in the ran- 
dom and scale free networks have the same second order 
degree. For values of t larger than r, the second order 
degree becomes larger for the random network than the 
scale free counterpart. This indicates that the second or- 
der degree in a scale free network tends to grow slower 
than in a random network. The log-log diagram in the 
inset, shows r(tj) in terms of tj for connection proba- 
bilities p = 0.9 and p — 0.6, which yields power laws 
r(tj) « tj a . 



It follows that 

A u = F u + G u A u -i 
In addition, we have to enforce that 

Ax=Fi, n k (i )Q = 
Developing the recurrence, we obtain 



(27) 



(28) 



A u — F u + G u (F u -i + G u -\A u -2) 

A u — F u + G u F u -i + G u G u -\{F u -2 + Gu — 2A U ^)... 



3. Second Order Degree Distribution 

We now consider the number of vertices with first order 
degree equal to k^ and second order degree equal to k^ 1 , 
which will be represented as N k (i) k (2). The obtention of 



, u u 



G u +i 



/=1 9=f 



(29) 



Expressed in terms of the original variables, it follows 
that 
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«fe(l)fc(2) — 



/2fc( 1 )+fc( 2 )+3\ I 



E 

/=i 



(fcW-iK (1) . 



if 



2fcW+/ + 2 



r(fc( 2 > + i)^ 



r(fe( 2 ) + 5) 



E/(/ + 1 )(/ + 2 )(/+ 3 w 



/=i 



1\ V 2fc (1) + ff + 3 



r 



By Eq.(21) we know the value of n/ (rif = Nf(t)/t). Substituting in the above expression: 
I 



2fcW+fc( 2 )+ 3\ x 



(2) r 

E 



2fc« + / + 2 



When fcW = 1, we have that 

2fc( 2 )(fc( 2 ) +7) 



g=f 



(fc(2) + l)(fc(2) + 2)(fc(2) + 3)(fc(2) + 4) 



(31) 



n ( 



2yfc(2)(fc(2) +7 ) 



2kW+g + 3j (fc(2) + l)(fc(2) +2)(fc(2) +3)(JfeP9 +4 



30) 



r 



This result is identical to the expression obtained by 
Krapivsky in for correlation degree. Developing the 
first part of Eq. (30): 



n k (i) k (2) 



n lk( 2) 



kW - 1 



2&(!) + fc(2) + 2 



r r(2fc(!) +/ + 2) 



E 



r 



r(fe(i) + /) 



"-fc(!)-l/ 



(32) 



Figure[fJ]shows the curves obtained for n fe <2) fe (2) for sev- 
eral values of feW and A:'- 2 -' varying from to 50. 

The total number of vertices with second order degree 
equal to fc^ 2 ) is given as: 



E n fe (i) fc <2) 

fc(l)=l 



(33) 



Figure[7|shows a log-log graph of the distribution iVwij 



(given by Eq.(21)) and iV^, (obtained from Eq.(33)). 
Note that the initial values in this distribution devi- 
ate from the linear relation typically observed for scale 
free networks. The second order distribution, given by 



.2) 



T- 



V. CONCLUDING REMARKS 

This work addressed the analytical characterization of 
hierarchical degrees of random and complex networks. In 



the case of random networks, we have shown that the hi- 
erarchical degree follows a power law, i.e. fc<»(f) ex t . 
By using the master equation approach, it has also been 
shown that the first order degree obeys a Poisson dis- 
tribution. Unlike the first order degree, the evolution 
of the second order hierarchical degree can not be easily 
described by a master equation as Eq.(6). This is a con- 
sequence of the fact that the adopted approach of unit- 
ing two connected components from the network into a 
single component implies several combinations of effects 
and respective terms to be incorporated into the master 
equation. In addition, we verified that the second order 
degree tends to grow slower in scale free networks than 
in random networks, with the position where these two 
values become equal following a power law. 

In the case of scale free models, namely the Barabasi- 
Albert network, it has been verified that the second or- 
der hierarchical degree, for a particular node, can be ex- 
pressed in terms of a logarithmic correction of the first 
order degree. Exact results have been obtained for the 
second order hierarchical degree. We observe that the 
generalization of such an approach to higher hierarchical 
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FIG. 6: Values of n fcC i )fc ( 2 ) for fc (1) = [1,5, 10, 15, 20, 25, 30, 35, 40, 45] and fc (2) in the range 0,50. 
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FIG. 7: Log-log distribution of P^ 1 ' and P^ 2 ' for Barabasi- Albert model of Scale Free Network. 



levels becomes substantially more complex because of the 
need to consider all recursions up to level h — 1. 

Possible continuations of the reported developments 
include the extension of the analytical expressions of hi- 
erarchical node degree to higher levels, as well as the 
derivation of analytical expression for other hierarchical 



features such as the hierarchical clustering coefficient and 
hierarchical number of nodes Q . 
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